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ABSTRACT 

The driving mechanism of the AGB superwind has become controversial in recent years. The 
efficacy of dust-driven mass loss has been queried. Spitzer observation of AGB stars in Local 
Group Galaxies show the surprising result that at low metallicity, AGB mass loss occurs at 
low luminosity, possibly lower than in the Galaxy, but only for carbon-rich stars. Oxygen-rich 
stars in the Galaxy and in lower metallicity galaxies have similar mass-loss rates only at high 
luminosities. To explain this dichotomy, we propose that the superwind has a dual trigger. 
The superwind starts either when sufficient excess carbon builds up for efficient formation of 
carbonaceous dust (which we propose occurs when Xqo — (C — 0)/Oq = 0.1), or when 
the luminosity reaches a value sufficient for a silicate-dust-driven wind (proposed at L = 
10 4 Z~ 4 / 3 L ). We show that this dual trigger fits the current observational constraints: the 
luminosity at which the superwind begins, and the predominance of carbon superwind star at 
low metallicity. We use stellar evolution models to check the consistency of our explanations 
and present detailed predictions of the luminosities at which the superwind is triggered for 
different metallicities and initial stellar masses. 
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1 INTRODUCTION 

The late phases of stellar evolution are crucial to galactic evolution. 
Stars in this phase eject much of their gas into space, enriched by 
the products of nuclear burning and newly formed dust. Of all the 
stars which have died in the Universe so far, over 90% are low and 
intermediate-mass stars (M ~ 1-8 Mq : LIMS). In consequence, 
LIMS are important contributors to the enrichment of the Interstel- 
lar Medium (ISM) by sources internal to the galaxy. 

LIMS eject between 20 and 80 per cent of their mass dur- 
ing the so-called superwind phase which terminates the high lumi- 
nosity asymptotic giant branch (AGB). This near-catastrophic but 
non-explosive wind removes effectively the entire hydrogen-rich 
envelope. During the superwind, dust forms in the ejecta, and the 
star tends to become self-obscured. After this final burst of activity, 
the star remains as a hot, compact white dwarf. The expanding cir- 
cumstellar ejecta become transparent, and briefly form an ionised 
planetary nebula, before dispersing into the interstellar medium. 

At the low effective temperature of AGB stars, the ejecta are 
largely molecular. The first stable and abundant molecule to form 
is CO: this locks up the least abundant of the C and O atoms. Other 
molecules and dust form from the remaining atoms. Dredge-up of 
primary carbon (produced by triple-a burning) causes the C/O ratio 
to increase over time. The chemistry changes dramatically when the 
C/O ratio passes unity, and this affects both the gas-phase species 
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(Millar et al. 2000, Willacy et al. 1997) and the dust composition 
(Treffers et al. 1974). In carbon stars, the C/O ratio (by number) 
exceeds unity, and carbon-rich molecules (e.g. C2H2, HCN) and 
carbonaceous dust result, whilst in oxygen stars (C/0< 1), metal 
oxides and silicate dust form. This dichotomy is not complete: in 
the inner shell of AGB stars shocks can lead to the formation of 
some molecules of the 'wrong' type (Decin et al. 2008). 

The mass-loss mechanism of AGB stars has become contro- 
versial. A two-step process is likely involved. First, shocks due to 
pulsations from the star extend the atmosphere (Jones et al., 1981) 
leading to dust formation. Second, radiation pressure on the dust 
drives an outflow; the gas is carried along through friction with 
the dust particles. Pulsations can drive mass-loss rates up to about 
lO _6 M0yr _1 ; the much higher rates observed are explained by 
dust-driven winds (Bowen & Wilson 1991). This simple distinc- 
tion between weak pulsation winds and dusty superwinds (Win- 
ters et al. 2000) has recently been disputed. Woitke (2006) find 
in their models that the dust opacity of silicates is insufficient to 
drive a wind, while for carbon dust the winds become very chaotic 
(Woitke & Niccolini 2005). For low-metallicity stars, dust-driven 
winds may fail due to the lack of dust (Zijlstra 2004): theoretical 
work by Bowen & Willson (1991) predicts that for metallicities be- 
low [Fe/H] = — 1 the wind becomes pulsation-driven. Mattson et al. 
(2008) argue that pulsation can drive a strong wind at low metallic- 
ity, set by the pulsation energy of the star. Recent theoretical models 
produce similar mass-loss rates for AGB stars with SMC, LMC and 
solar metallicities (Wachter et al. 2008). 
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The effect of dust should be most noticeable when compar- 
ing stars at different metallicity. This has become possible with 
the Spitzer Space Telescope (Werner et al. 2004), which provided 
the sensitivity necessary to obtain mid-infrared spectra and derive 
mass-loss rates for the whole range of AGB masses and luminosi- 
ties in the Magellanic Clouds and a few other Local Group Galax- 
ies. These galaxies trace in general metal-poor stars with a range of 
metallicity and ages, and the known distances means that absolute 
mass-loss rates and luminosities can be established. The expecta- 
tion for reduced mass loss efficiency at low Z has not been con- 
firmed. Instead, the Spitzer surveys have shown that high mass loss 
rates with dusty winds can be reached even at very low Z, as seen 
in the Fornax dSph galaxy at [Fe/H]~ —1 (Lagadec et al. 2008). 

One unexpected result of the Spitzer surveys is the almost 
complete dominance of carbon stars among the mass-losing stars 
in these low metallicity galaxies. The estimated mass-loss rates for 
these carbon stars are very similar to the ones observed in Galac- 
tic AGB stars, while the few oxygen-rich AGB stars observed have 
mass-loss rates one or two orders of magnitude lower than their 
Galactic counterparts (Sloan et al. 2008). All this seems to indicate 
that the abundance of carbon and oxygen is an important parameter 
for the mass-loss from AGB stars at low metallicity. We will dis- 
cuss in the present letter at what point during the stellar evolution 
the superwind begins, and how it is triggered. 



2 THE AGB STARS SAMPLE, LUMINOSITY AND 
MASS-LOSS RATES 

2.1 Sample selection 

The stars from our sample were selected as reddened AGB stars in 
a 2MASS Mfcvs J-K diagram. FigQ] shows the colour-magnitude 
diagrams for the Large Magellanic Cloud, Small Magellanic Cloud 
and Sagittarius dwarf spheroidal galaxy (hereinafter LMC, SMC 
and Sgr dSph). For the LMC and SMC, the dots indicate optically 
identified carbon stars. For Sgr dSph, no such catalogue is known. 
Instead we selected all stars within a 1 -degree region towards the 
galaxy, and removed any star with a close match in (K, J— K) with 
a neighbouring field. We finally used a colour criterion to remove 
the remaining sequences corresponding to foreground stars, to the 
blue of the Sgr sequence. The majority of the remaining stars will 
be oxygen-rich RGB stars: the AGB is visible as the less populated 
upper range of the sequence. In all three plots, open circles indi- 
cate stars selected for Spitzer observations. For Sgr dSph, these are 
selected from a larger area than the plotted unreddened stars. 

All three systems show a range of metallicities. We take as 
representative values [Fe/H] = —0.3 for the LMC, —0.7 for the 
SMC and —0.55 for Sgr dSph. The distance moduli are taken as 
18.54, 18.93 and 17.02, respectively. 

Several different Spitzer programs have obtained mid-infrared 
spectroscopy of AGB stars in the Local Group. The description 
below refers mainly to works published in Zijlstra et al. (2006); 
Matsuura et al. (2006); Groenewegen et al. (2007); Lagadec et al. 
(2007) and Sloan et al. (2006; 2008) for carbon rich stars in the the 
Magellanic Clouds (27 stars in the LMC and 33 in the SMC) and 
van Loon et al. (1999) for the oxygen-rich stars in the LMC (26 ob- 
jects). Matsuura et al. (2007) for 3 AGB stars Fornax and Lagadec 
et al. (2008) for 6 others in the Sagittarius Dwarf spheroidal galaxy. 
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Figure 1. J— K vs Mk for the LMC and SMC samples. The dots indicate 
2MASS field stars in each galaxy. For the LMC and SMC, all field carbon 
stars are shown. For Sgr dSph, a cleaned field star population (oxygen rich 
or carbon rich) towards the central square degree of the galaxy are shown: 
see text. 



2.2 Luminosities 

The important onset of the superwind is best traced by the luminos- 
ity distribution of the dusty stars, compared to those of unreddened 
AGB stars. 

Fig. [2] shows the distribution of the AGB carbon star popula- 
tions in the LMC and SMC (bottom panels; left and middle). These 
were taken from spectroscopically identified catalogues of carbon 
stars. The bolometric magnitudes were derived using the J — K, K 
magnitudes: the bolometric correction for carbon stars was taken 
from Whitelock et al. (2006). It is important to note the uncer- 
tainties inherent in such bolometric corrections. Most of the near- 
infrared spectral region where the stellar energy distribution peaks 
is hidden by our atmosphere. In addition deep molecular bands are 
exhibited by the stars. Bolometric corrections derived from oxygen- 
rich stars would in general overestimate the luminosity of carbon 
stars, as the former corrects for deep water bands which coincide 
with opaque spectral regions in our atmosphere. 

The histograms of the absolute bolometric magnitudes distri- 
butions for the dusty stars observed with Spitzer, are shown in the 
upper panels of Fig. [2] For these, the luminosity was obtained by 
integrating over the Spitzer spectra and broadband colours. 

We also show some corresponding data for Sgr dSph. The 
panel shows stars selected from the sample of Fig. [I] using a colour 
criterion to select the reddest branch seen in the figure. This will be 
dominated by oxygen-rich stars. The observed distribution clearly 
shows the dominant inclusion of RGB stars. We also do not know 
how many of the AGB stars are carbon-rich. A comprehensive sur- 
vey for carbon stars in Sgr dSph is clearly warranted. 

Table Q] gives some representative parameters for the carbon 
stars in each galaxy, derived from the distributions in Fig. [2] The 
second column lists the magnitude where the luminosity function 
of the unreddened carbon stars becomes flat: this identifies the 
point where the large majority of AGB stars has become carbon- 
rich. The third column indicates the life time of the carbon star 
phase: this is derived from the width of the flat part of the lu- 
minosity function, converted to a time scale using the relation 
dM hol /dt = -8.25 x 10~ 7 mag yr" 1 (Wood 1990). The last col- 
umn indicates the bolometric magnitude at which the onset of the 
superwind occurs: this is taken as the point where the luminosity 
distribution of the mass-losing stars reaches its maximum. These 
numbers should be used with caution: especially the last column is 
affected by low-number statistics. 

For the LMC, the table and figure show that the mass loss 
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Table 1. Approximate numbers for life times and magnitudes of carbon 
stars, derived from the current samples. Afl" 1 / indicates the magnitude 
where the large majority of stars have become carbon-rich: i.e, the luminos- 
ity function becomes flat. M bo i(gw) indicates the onset of the superwind. 
Life times of the carbon-rich phase are taken from Lagadec et al. (2007). 
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Figure 2. Bolometric magnitudes for the observed stars (top panels) and for 
field stars. For the LMC and SMC, only field carbon stars are included. For 
Sgr dSph, the field population is for the red colour sequence only, which 
includes oxygen-rich stars. 



starts shortly after the star becomes carbon rich. For the SMC, there 
is a longer delay. Also, for the LMC, most stars become carbon-rich 
late on the AGB. For the SMC, the shift in chemical type coincides 
with the tip of the RGB, which is approximately the same as the 
onset luminosity for thermal pulses. Thus, the SMC stars become 
carbon-rich very soon after the thermal pulses begin. 

The last column, surprisingly, seems to indicate that the mass 
loss starts at a lower luminosity, i.e. earlier on the AGB, for the 
SMC than it does for the LMC and the Galaxy (M bo i ~ -5.5). 
One should note that the samples were not selected in quite the 
same way, and so any conclusion would require more work. 



2.3 Mass-loss rates 

The mass-loss rates of stars observed with Spitzer have been esti- 
mated by fitting spectra using a 1-d radiative transfer code, includ- 
ing amorphous carbon dust and silicate carbide dust. The fitting 
process and results are described in Groenewegen et al. (2007) for 
the MCs stars, Matsuura et al. (2007) for the Fornax stars and La- 
gadec et al. (2008) for Sgr dSph. 

The results for the LMC, SMC, Sgr dSph and Fornax are 
shown in Fig(3] We added the O-rich stars sample from van Loon 
et al. (1999) and the well-studied Galactic carbon star IRC+10216. 
The figure shows that the mass-loss rates are highest for the LMC 
stars, and tend to be lower for the SMC. This is under the assump- 
tion of a constant gas to dust ratio. There is a tentative indication 
that the gas-to-dust ratio is higher at the lowest metallicities (Mat- 
suura et al. 2007). 

The mass-loss rates are around 10 -5 Mq yr -1 which is in the 
same range as the high mass-loss stars in the Galaxy. A few Galac- 
tic stars may reach rates up to ten times higher, but even in the 
Galaxy, such stars are very rare. Overall, there is no evidence that 
peak (gas) mass-loss rates depend on Z for carbon-rich stars. 



2.4 Carbon star formation 

One unexpected result of the Spitzer surveys is the almost com- 
plete dominance of carbon stars among the mass-losing stars in the 
Magellanic Clouds. The original selection criteria did not separate 
the two classes, so that this result shows a real effect. Part of the 
explanation is that, at lower metallicity, more stars become carbon 
stars. Less dredge-up is required to overcome the original oxygen 
abundance, and acquire the C/0> 1 needed to form a carbon star. 



Based on the luminosities, the progenitor masses of the mass-losing 
carbon stars are ~ 1.5-2.5 Mq (Zijlstraet al. 2006). 

The picture that emerges from this is that, at LMC metallic- 
ity, all stars in the ~ 1.5-2.5 Mq range are C stars by the time they 
develop substantial mass loss rates and they remain C stars until 
their AGB evolution is terminated. The more luminous and mas- 
sive, but rare, OH/IR stars in the LMC ([e.g. Marshall et al. 2004) 
can survive as O-rich objects well into the dust-enshrouded super- 
wind phase, possibly to the end of their AGB evolution. 

Table Q] lists the bolometric magnitude at which the onset of 
the superwind occurs. This happens at a somewhat lower luminos- 
ity (L ~ 5 10 3 Lq) than is found in our Galaxy (L ~ 10 4 Lq; 
Afbol ~ —5.5). The difference between the LMC and the Galaxy 
is that the former become carbon stars at a point on the AGB where 
the latter are still oxygen-rich. We interpret this as an indication 
that carbon-rich stars are more efficient at mass loss, even where 
the oxygen-rich stars are more metal-rich. The finding agrees with 
(but is less extreme than) the prediction of Woitke (2006). 

The few oxygen-rich stars in our MCs sample show mass-loss 
rates 1-2 orders of magnitude less than the carbon stars (Sloan et 
al., 2008). Thus, the mass-loss efficiency in oxygen-rich stars is 
more affected by metallicity. This may not be unexpected: oxygen- 
rich dust depends on metallicity-limited elements (Si, Al), while 
amorphous carbon depends on self-produced carbon, and therefore 
the dust abundance formation for oxygen-rich stars is much more 
strongly dependent on metallicity. Taking the argument in reverse, 
the implication would be that dust plays a role in driving the super- 
wind of oxygen-rich stars in the Galaxy. 

To estimate the evolution of the carbon and oxygen abun- 
dances in the circumstellar envelopes of the AGB stars we ob- 
served, we used the stellar structure models described in Karakas 
& Lattanzio (2007) and made available for download. It calculates 
yields from detailed models of low and intermediate-mass AGB 
stars with masses in the range 1-6 Mq . Models for four different 
metallicities are available: Z = 0.02 (solar), Z = 0.008 (LMC), 
Z = 0.004 (SMC) and Z = 0.0001. We have thus studied the evo- 
lution of the C/O abundance ratio with the luminosity of the star 
for these four metallicities. At lower metallicity, the initial oxygen 
abundance is lower and the stars become carbon-rich earlier in their 
evolution. At solar Z, stars with initial masses below 3 Mq never 
become carbon-rich, while in the LMC and SMC, this limit drops to 
1.75 and 1.5 Mq respectively. For a star with a 2.5 Mq initial mass, 
the C/O ratio becomes larger than 1 at L ~ 9.5 10 3 Lq in the LMC, 
~ 6.5 10 3 L Q in the SMC and ~ 7.5 10 3 L Q for Z = 0.0001. (In 
the last case, the larger value is due to a metallicity effect in the 
core-mass-luminosity relation). The evolution of the C/O ratio with 
luminosity for the four different metallicities are shown in Fig[3] 
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Figure 3. Mass-loss rates as a function of luminosity for different Local Group AGB stars. Panels a,b,c and d show the evolution of the C/O ratio with 
luminosity using the models described in Karakas & Lattanzio (2007.). Values of the different initial masses in solar units are on top of each panel. 



At low metallicity, third dredge-up is more efficient in making 
stars carbon rich: the same amount of primary carbon will have a 
larger effect on more metal-poor stars. This process does not cease 
once the star has become carbon rich: further dredge-up will con- 
tinue to enhance the C/O ratio. We may therefore expect that metal- 
poor carbon stars have a higher C/O ratio than do metal-rich carbon 
stars. Matsuura et al. (2005a) find evidence that LMC stars have 
C/O ratios of ~ 1.5, versus ~ 1.1 for typical Galactic carbon stars. 
Ratios in lower metallicity systems (e.g. SMC) could be higher. 

This immediately affects the potential amount of amorphous 
carbon dust, which will depend on the amount of 'free' carbon. 
Instead of the C/O ratio, a more relevant parameter is 

C-O 



X c = 



o. 



which measures the excess carbon by number, scaled to the solar 
oxygen abundance. Xc rapidly increases towards lower metallicity. 
The corresponding increase in the C2H2 abundance is in fact seen, 
with the acetylene bands becoming much stronger (larger equiv- 
alent width) in the sequence Galaxy-LMC-SMC (Lagadec et al. 
2007). In contrast, the absence of HCN bands suggests its abun- 
dance becomes limited by nitrogen, which becomes enhanced only 
in the more massive stars. 

The dust formation efficiency may depend primarily on Xc- 
(The seeds required for the dust formation (e.g. TiC) may intro- 
duce a dependency on metallicity.) For mass loss driven by radia- 
tion pressure on carbonaceous dust, we make the assumption that 
the superwind is triggered when a limiting value for Xc is reached. 
For Galactic stars, this occurs for C/0= 1.1 (Xc = 0.1) (Ohnaka 
et al., 2000). We predict that for the LMC the superwind begins at 
C/0= 1.25 and for the SMC, at C/0= 1.5. The stronger acetylene 
bands at lower metallicity, and the indicative C/O ratios (Matsuura 
et al. 2005a) suggest similar but possibly slightly higher C/O ratios. 



3 DISCUSSION 

For oxygen-rich stars, Galactic stars show a superwind for lumi- 
nosities of ~ 1O 4 L0. The mass-loss rate of a dust-driven wind 
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Figure 4. Left:Initial mass of the stars as a function of the luminosity at 
which their C/O ratio becomes larger than 1 . Right: Distribution of the ini- 
tial mass of the AGB stars in our LMC (black line)and SMC sample (dashed 
line) as derived from the relation son the left. Values are in solar units. 

scales as L B > 4 (van Loon, 2007) and is proportional to the dust-to- 
gas mass ratio *]/. As silicate dust depends on Si and O, which are 
not enriched in AGB stars, we assume that "if scales as Z. We now 
make the assumption that the superwind starts when this relation 
predict a critical mass-loss rate. Scaling from the Galaxy, we then 
expect that an oxygen-rich dust-driven superwind (OSW)(larger 
than lO _7 M0yr _1 ) is triggered for 

-4/3 
L 



Losw — 10 ( -=- 

For the LMC and SMC this corresponds to 3.4 x 10 4 and 8.5 x 
10 L©, with an uncertainty of a factor of ~2 due to the estimation 
of the value of 10 4 L© and the metallicity dispersion within the 
galaxies . Fig.[3]shows that the least luminous OSW star in the LMC 
has L ~ 1.8 10 Lq. At lower luminosities, oxygen-rich stars have 
much lower mass-loss rates than do carbon rich stars. 

For luminosities below this limit, a superwind is seen only for 
carbon-rich stars (Fig[3j. There is no clear evidence of increasing 
mass-loss rates with increasing luminosity. The carbon-rich super- 
wind (CSW) in low Z AGB stars occurs even at luminosities below 
1 x 1O 4 L0, showing the efficacy of carbon dust in driving mass loss. 

We thus propose that the superwind is triggered when either 
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Table 2. The luminosities for which C/O reaches 1, or Xc 0.1 for different metallicities and initial masses, derived from the stellar evolution models, 
indicates that no model was available and — indicates that the considered star never become carbon-rich. Luminosity are in unit of 1O 4 L0. 
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i) the luminosity of an AGB star becomes larger than a lumi- 
nosity limit (Lii m it ~ W 4 (Z/Zq)~ 4 ' 3 L©), regardless of chem- 
istry, over which the radiation pressure becomes dynamically im- 
portant. This is the dominant mechanism at solar metallicity 

ii) or the amount of free carbon reaches Xc ~0.1, when a 
CSW starts. This is the dominant mechanism at subsolar metallic- 
ity, for stellar masses above the limit for third dredge up. 

Pulsation-driven winds are not included here: we expect these 
may become important for low-metallicity, low-mass stars which 
do not experience third dredge-up. 

Fig|4] shows the relation between the initial mass of the star 
and the luminosity at which Xc ~0. 1 according to the stellar mod- 
els. If the superwind is triggered when Xc ~0.1, then we can esti- 
mate the initial mass of the observed stars using this relation. The 
histograms of the observed stars in the LMC and SMC are shown in 
FigS These histograms show that most of the observed stars have 
an initial mass between 1.5 and 3.5 Mq. If Xc is the main factor 
to trigger the superwind, then we can also predict the luminosities 
at which the superwind will be triggered for different masses and 
metallicities. These values are shown in Table|2] 



4 CONCLUSIONS 

In this letter, we made use of a Spitzer spectroscopic survey of 
mass-losing AGB star in Local Group low metallicity galaxies to 
study the superwind from the stars. We have shown that the lumi- 
nosity at which the superwind is observed is decreasing with lower 
metallicity for carbon-rich stars. Oxygen-rich stars in the Galaxy 
and in lower metallicity galaxies have similar mass-loss rates for 
similar luminosities. To explain this dichotomy, we propose that the 
main agent to trigger the superwind in low metallicity carbon stars 
is the abundance of carbon. When the C/O ratio becomes larger 
than unity, carbon-rich dust forms and the radiation pressure on 
this dust becomes more efficient. We used stellar evolution models 
to check the consistency of our explanations and made some predic- 
tion of the luminosities at which the superwind should be triggered 
for different metallicities and stars initial masses. 
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